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ABSTRACT 

We present magnet ohydro dynamic numerical simulations of the late post-supernova hypercritical 
accretion to understand its effect on the magnetic field of the new-born neutron star. We consider as 
an example the case of a magnetic field loop protruding from the star's surface. The accreting matter 
is assumed to be non magnetized and, due to the high accretion rate, matter pressure dominates over 
magnetic pressure. We find that an accretion envelope develops very rapidly and once it becomes 
convectively stable the magnetic field is easily buried and pushed into the newly forming neutron 
star crust. However, for low enough accretion rates the accretion envelope remains convective for an 
extended period of time and only partial submergence of the magnetic field occurs due to a residual 
field that is maintained at the interface between the forming crust and the convective envelope. In 
this latter case, the outcome should be a weakly magnetized neutron star with a likely complicated 
field geometry. In our simulations we find the transition from total to partial submergence to occur 
around M ~ IOMq yr~^. Back-diffusion of the submerged magnetic field toward the surface, and the 
resulting growth of the dipolar component, may result in a delayed switch-on of a pulsar on time-scales 
of centuries to millenia. 

Subject headings: accretion — magnetohydrodynamics (MHD) — magnetic fields — stars: neutron 
— supernovae: individual (SN 1987A) 



1. INTRODUCTION 

Most observed neutron stars show clear evidence for 
the presenc e of strong magnetic fields. In the case of 



magnetars (Woo ds fc Thompson] 2006) the estimated 
strength of the surface magnetic field is of the order of 
10^^ G while for the grind of the mill radio pulsar 10^^ 
G is a typical value. Lower magnetic fields are, however . 



found in millisecond pulsars (iPhinney fc Kulkarni] 1994t 

m 



and i n neutron stars in low-mass X-ray binaries ( [Fsaltis 
|2006[ ), but in both cases past or present, respectively, ac- 
cretion is thought to be the cause of the magnetic field 
reduction. There is, however, a small group of neutron 
stars, found in young supernova re mnants and dubbed 
CCOs ("Central Compact Object" , |Pavlov et al.|[2QQ2t , 
which exhibit little or no evidence tor the presence of a 
magnetic field. 

The case of the supernova SN1987A is also intriguing 
since no evidence for the presence of a compact object 
has yet been foun d, despite extensive searches (see, e.g.. 



Manchester 



120071). The presence of a young energetic 
pulsar is clearly excluded, but a slowly rotating neutron 
star, with period P > 100 ms, and with a not too strong 
dipolar magnetic field, B d < 10^^ G, is still compatible 



with all current dat a (Manchester||2007[ |Ng et al.||2009 



Larsson et al. 2011). Such a period and low magnetic 



field are within the range of values inferred, when p os- 
sible, in some of the CCOs ( [Gotthelf fc Halpern[[2Q08) . 
The possibility that this supernov a may have proauced a 
black- hole, as proposed, e.g., by Brown fc Bethe |l994| 

oF 31 



on the basis of a very soft dense matter equation of state 
and, consequently, a low neutron star maximum mass 
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around I.GMq, is now very slim i n view of the existence 
of a 2Mq pulsar (Demorest et aTj lgOlOl). 
The origin of neutron star magnetic fiel ds is still an 



unsol ved problem (for recent reviews, see jSpruit 2008 
2009). Two main mechanisms, a fossil field from the pro 



genitor compressed during the core-collapse and a proto- 
neutron star dynamo, are still competing and are likely 
both needed to explain the large variety of observed field 
strengths. In these scenarios, the magnetic field gen- 
eration and/or adjustment process terminates within a 
minute after the neutron star's birth. 

After this early field development the story is not nec- 
essarily over. The supernova shock is still pushing its 
way through the outer layers of the progenitor and, if it 
encounters a density discontinuity, a reverse shock can 
be generated. Depending on its strength and on how 
far out it was generated, this reverse shock can induce 
strong accretion onto the new-born neutron star on time 
scales of hours. Notice that this delayed accretion has to 
be distinguished from the initial fall-back which occurs 
seconds after the initial core collapse. A particularly fa- 
vorable configuration for such late accretion is present in 
the cases the progenitor had a tenuous H/He envelope 
surrounding a de nse He core, as in Type lib supernovae 
or in SN 1987A ([S nia£tt][2Q^ . 

Chevalier ( 1980p argued in favor of such late accretion 
the case or SN 1 9 8 7A and developed a simple a nalyt i- 



cal model for it (see also Brown fc Weinga rtner^994 for a 



similar model and Figure ernal et al. 2010 for a de 

piction of these scenarios). Highly super-Eddington ac- 
cretion results, in which the photons are trapped within 
the accretion fiow and the energy liberated by the accre- 
tion is lost through neutrino emission close to the neutron 
star surface. Such a regime has been termed "hypercrit- 



ical accretion" ([Blondiijljj^ |Chevalier|p89l |Houck 



Chevalier||1991[) and requires an accretion rate rh higher 
than about 10 rhEdd, where rfiEdd ^ 10^ g cm~^ s~^ is 
the Eddington rate. We note that this regime is rel- 
evant n ot only for SNe, but for Gamma Ray Bursts 
as weh ([Narayan fc Quataert | [2QQ5l | Nakar||2QQ7[ |Lee fc 
Ramirez-^uiz 2007 ; Gehrels et alj 2009 ), allowing rapid 

±1 



mass accretion onto newborn black holes to produce the 
required power to account for isotropic equivalent lumi- 
nosities that can reach a bove 10^^ erg s"-*^ in the prompt 
emission (Narayan et aL||2001 ; Kohri fc Mineshige 2002; 
' Lee et aL 2005|), and possibly extended emission episodes 
as we ll ( [kosswog 2007 Lee et al. 2009 Metzger et al. 



2010). In the case ot core collapse events, after the re- 



verse shock hits the neutron star surface a third shock 
develops and starts moving outward against the infalling 
matter. Once this accretion shock stabilizes it will sep- 
arate the infalling matter from an extended envelope in 

quasi-hydrostatic equilibrium.^ 

Following the suggestion of [Muslimov fc Page ( 1995"), 
Geppert et al. (11999) presented simple ID ideal JVLHD 
simulations of the effect of this post-supernova hyper- 
critical accretion on the new-born neutron star magnetic 
field. The result was a rapid submergence of the field 
into the neutron star. After accretion stopped, the field 
could diffuse back to t he surface and result in a delayed 
switch-on of a pulsar (Michel 1994 Muslimov fc Page 
1995 ). Depending on the aniount oi accreted matter, 
the submergence could be so deep that the neutron star 
may appear a nd remain unmagnet ized for more than a 



Hubble time ([Geppert et al. 



recently revisited by'Ho"(201 



field evolution of the CCOsT 



1999). This scenario was 
and applied to study the 



(1999) 



The ID MHD simulations of Geppert et al. 
were, however, carried out wit h some considerab le sim- 
plifications. In a previous work ( Bernal et aL|201Q , Paper 



I hereafter) we presented results of 2D MHD simulations 
of hypercritical accretion onto a magnetized neutron star 
surface, with simple magnetic field configurations. We 
simulated an hypercritical accretion flow in a rectangular 
domain, with the neutron star surface lying at its base, 
in which an initial uniform, horizontal or vertical, mag- 
netic field is present, as well as an horizontal field with 
strength decreasing with height. For highly- hypercritical 

accretion rates, rfi ^ lO^^mEdd, complete submergence 
of the magnetic field was observed in a time scale of a 
few tens to a few hundreds of milliseconds. 

In the present work we extend the scenarios considered 
in Paper I, to include an initial magnetic field configura- 
tion which consists of a magnetic loop protruding from 
the neutron star surface. We mostly consider the 2D 
case in a cylindrical box and one example of a 3D model. 
We also extend our simulations to lower accretion rates, 
down to lO^mEdd- With the field strengths we employ, 
matter and ram pressures still strongly dominate over the 
magnetic pressure. Moreover, we assume the accreting 
matter is unmagnetized. 

In ^we describe the numerical method and a simple 
analytical model. In ^we describe our results, of which 
an interpretation is proposed in Q Finally, we conclude 
in© 



2. THE NUMERICAL HYPERCRITICAL 
ACCRETION MODEL 

The present work is a direct extension of our previ- 
ous work presented in Paper I. As described in Paper I, 
we use a customiz ed version of the hy drodynamic code 
AMR FLASH2.5 ( [Fryxeh et an[2QQQf , with the Split 8 
wave solver which solve the whole set of MHD equations. 
We work in the ideal MHD regime with only numerical 
resistivity and viscosity. The matter equation of state is 
an adaptation of FLASH'S HELMHOLTZ package which 
includes contributions from the nuclei, e~ — e+ pairs, and 
radiation, plus the Coulomb correction. Neutrino energy 
losses are dominated by the e~ — e+ annihilation process, 
but we also include the photo- neutrino, plasmon decay, 
e-ion bremsstrahlung and synchrotron, processes, which 
are implemented in a customized module (see Paper I). 
No nuclear reactions are taken into account. 

Since the MHD equations can be solved by FLASH 
only in cartesian coordinates, we consider wide columns 
in 2D and 3D, with a base of Ax = 2 • 10^ cm centered 
on X = in 2D, or Ax x Az = {2 - 10^) x (2 • 10^) 
cm^ centered on (x, = (0,0) in 3D, and a height Ay 
of several times 10^ cm, with y = 10^ being the neutron 
star surface. The mapping of the surrounding cone above 
the neutron star into a column is illustrated in Figure [l] 
We consider as magnetic initial condition a magnetic field 
loop, in the shape of an hemi-torus. On the central hemi- 
circle of the loop the field has a strength Bq = 10^^ G 
and about it is shaped as a Gaussian, i.e., with strength 
B{d) = B^exp [—{d/RL)'^)]^ d being the distance to the 
loop central hemi-circle and Rl = I km. The two feet 
of the loop are centered at x = — 5 and x = +5 km, and 
z = in the 3D model. The gravitational acceleration is 
taken as gy = —GM/y'^ and we assume a neutron star 
mass of IAMq. 

As boundary conditions, we impose mass inflow along 
the top edge of the computational domain, and periodic 
conditions along the sides. At the bottom, on the neu- 
tron star surface, we use a custom boundary condition 
which enforces hydrostatic equilibrium (see Paper I). For 
the magnetic field, the two lateral sides are also treated 
as periodic boundaries, while at the bottom the field is 
frozen from the initial condition, i.e., the two feet of the 
loop are anchored into the neutron star and no field can 
be pushed into the star by the accretion. On the top 
boundary the magnetic field is set to zero, i.e., we as- 
sume the accreting matter to be non-magnetized. 

We are interested in following the evolution of the mag- 
netic field under the heavy hypercritical accretion from a 
reverse shock (see Figure 1 in Paper I) which reaches the 
neutron star surface hours after the core collapse. The 
various accretion rates, m, we consider will be expressed 
in terms of a fiducial rate 



mo 



1.75 X 10^^ gcm" 



^s-i 



(1) 



and always assumed to be constant during the entire sim- 
ulations. This corresponds to a total accretion rate onto 
the neutron star of 



Mo 350 Mq yr" 



(2) 



and is the accretion rate originally estimated by |Cheva^ 
lierl (11989') for SN 1987A. We consider two different sce- 



narios regarding the initial conditions for the accreting 
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matter. In the first case, corresponding to free fall, we 
start the simulation just before the reverse shock reaches 
the magnetic loop, which is initially immersed in a low 
density, 10^ g cm~^, medium, and above it the falling- 
back matter has a free-fall velocity and a density ob- 
tained by mass conservation: 



Vf[ — \/2GM/y and = rn/vf[ 



(3) 



As we will see, after a few tens of milliseconds an ac- 
cretion envelope develops which is in a quasi-hydrostatic 
equilibrium, separated from the continuous accretion in- 
flow by an accretion shock. 

In the second set of scenarios, we start with this quasi- 
hydrostatic equilibrium envelope (QHEE hereafter) and 
the magnetic loop immersed within it. The structure of 
the QHEE can easily be obtained analytically, 

P = Psh{ysh/yf, p = pshivsh/yf^v = Vsh(^/^sh)^, (4) 

where the first two values come from imposing hydro- 
static equilibrium for a poly tropic equation of state, 
P (X p^ with 7 = 4/3, and the velocity is fixed by mass 
conservation. Once the shock position y^h is known (see 
below for its determination), Pgh and psh are determined 
by the strong shock condition and v^h by mass conserva- 
tion as 

72 I 

Psh = ypsh^sh^ Psh = 7pff(^sh), ^sh = ;^^ff(^sh). (5) 

For given M and and a fixed m, the vertical location of 
the accretion shock, ^gh, is controlled by energy balance 



GMrfi 
R 



POO 

= / ^u{y)dy 

JR 



(6) 



between the accretion power and the integrated neutrino 
losses, per unit neutron star surface area. Neutrino losses 
are dominated by e~ — pair annihilation for which a 



simple rate is ( Dicus|p^72 ) 



1.83 X 10 



-34 d2.25 



erg cm s 



(7) 



Due to the resulting strong y dependence of e^^, y~^jthe 
value of the upper limit in the integral of Equation (pi) is 
not important and this fixes the height of the accretion 
shock as 



ysh-Sx 10^ (mo/m) 



10/(7-9) 



cm . 



(8) 



A stationary envelope, in quasi- hydrostatic equilibrium, 
will expand, or shrink, so that the physical conditions 
at its base allow neutrinos to carry away all the en- 
ergy injected by the accretion. Once emitted, neutri- 
nos will act as an energy sink provided the material is 
optically thin to them. Under the present conditions 
of density and temperature at the base of the flow, the 
fluid consists mainly of free neutrons, protons and elec- 
trons. The main sources of neutrino opacity are then 
coherent scattering off neutrons and protons and pair 
annihilation. For example, t he corresponding cross sec- 
tion for coherent scat tering CTubbs & Schramm" 1975; 
Shapiro fc Teukolsky|l9 83) is = (l/4)c^o[^^/(meC^)]^ 

. As these are thermal 
/cbT, and with temper- 



where cTo = 1.76 X 10-^"^ c] 
neutrinos, their energy is Ej^ 



Accretion flow 




Fig. 1. — Mapping of a portion of a spherically symmetric accre- 
tion flow onto a neutron star into a 2D, or 3D, cartesian domain. 
Illustrated is a magnetic fleld loop that will react, anisotropically, 
to the accretion. 



^ 7 xlO cm^. The maximum densities reached at 



atures T^IO K^IO MeV as we will find, we have 



the bottom of the envelope will be below 10^^ g cm ^, 
and in such conditions, the neutrino mean free-path is 

ly = (^nctn)"^ ^ 2.5 X 10^ cm, which is safely larger 
than the depth of the dense envelope, of the order of a 
few km. Above this dense region the envelope density 
decreases rapidly. Equation Q , and the whole envelope 
is practically transparent to neutrinos. We will, hence, 
ignore neutrino absorption and heating. 

When using a QHEE as initial condition, the P, p, 
and V profile of Equation Q are modified in the first 
two kilometers just above the neutron star surface to 
smoothly turn into raising P and p, and decreasing v 
power-law profiles. 

Besides been used as an initial condition in the QHEE 
scenarios, this analytical envelope model will be useful in 
the free-fall scenarios for comparison with the numerical 
results when the system as reached a stationary state. 

In Paper I we performed detailed comparisons of the 
hydro solver and the MHD solver with zero magnetic 
field and obtained excellent agreement, as well as agree- 
ment with the above described analytical envelope mod- 
els when a stationary states had been reached. More 
details can be found in Paper 1. 

3. RESULTS 

We will here present results from two series of sim- 
ulations, with different accretion rates. In the first se- 
ries, we consider highly hypercritical accretion rates, 
m = (1, 10, 100) • mo, while in the second series weakly 
hypercritical rates, rh = (10~^, 10~^, 10~^, 10~^) -rho are 
studied. We will describe the scenario of matter initially 
in free-fall, with the resulting formation of an accretion 
shock and the development of a QHEE, and the simpli- 
fied scenario where we start with a previously formed 
QHEE and follow its evolution. With the high accretion 
rates of the first series we are able follow the complete 
evolution of the accretion shock and the development of 
the QHEE in the free-fall scenario and compare the re- 
sults with the simplified purely QHEE scenario. In con- 
tradistinction, in the presence of the lower accretion rates 
of the second series, the accretion shock would move out- 
ward to very large distances and leave the computational 
domain and we are, hence, forced to restrict ourselves 
to the initial QHEE scenario. Our distinction between 
highly and weakly hypercritical rates is thus purely nu- 
merical and does not implies different physical conditions 
and/or results. At the highest accretion rate, 100 mo, 
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since the accretion shock expands to smah heights we 
win present results of a 3D MHD simulation while for the 
other cases only 2D simulations have been performed. 

3.1. Highly hypercritical accretion rates 

We describe here results for highly hypercritical accre- 
tion rates, m = (1, 10, 100) • mo, for which we are able to 
follow the evolution of the accretion shock. In Figure [2] 
we plot their final radial profiles of density, pressure, ve- 
locity and temperature for 2D simulations in the free-fall 
scenario. Here the final state is defined once the accre- 
tion shock, and the included envelope, have reached a 
quasi-stationary state. In the density and pressure pro- 
files the piling up of matter close to the neutron star 
surface is notorious (an effect not accounted for in the 
analytical approach). Notice that, although the system 
has reached a quasi- stationary state in all cases, signifi- 
cant noise remains in the velocity profile: because of the 
periodic boundary condition on the vertical sides mat- 
ter can freely fiow in the horizontal direction, as well 
as bounce off the neutron star surface, which prevents a 
full stationary state from being reached. Nevertheless, 
the mean velocity profile is close to the analytical one, 
indicating that the system has largely relaxed despite 
these fiuctuations. In these free-fall simulations, for our 
fiducial accretion rate mo the system reaches the quasi- 
stationary state in about 600 ms, whereas for higher ac- 
cretion rates this time is substantially reduced: 300 ms 
for lOrho and 100 ms for 100 mo- 

We show in Figure [3} for a 100 mo accretion rate, color 
maps of density with magnetic field contours superim- 
posed, for the 2D free-fall scenario. The initial panel at 
t = shows the reverse shock just before it reaches the 
magnetic field loop. When the reverse shock hits the NS 
surface and bounces a violent, convective, layer appears, 
the accretion shock develops and rapidly moves upward 
against the infalling matter {t = 1 ms panel). The fiow 
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Fig. 2. — Radial profiles of density, pressure, velocity, and tem- 
perature, for highly hypercritical accretion rates m = (1, 10, 100) • 
mo (labeled), in the 2D free-fall case, once the stationary state has 
been reached. 



smoothens gradually as the accretion shock stabilizes. 
Instabilities of the Ray leigh- Taylor type are present in 
this regime, but they disappear when the system reaches 
equilibrium (see the panel at t = 100 ms). The mag- 
netic loop is immediately torn apart by the reverse shock 
and the initial strong convection which is dragging the 
field with it within the developing envelope (panel at 
t = 1 ms). Although it has a turbulent dynamics, the 
magnetic field, having its feet frozen into the neutron 
star surface, always remains confined below the accre- 
tion shock. After 50 ms, the magnetic field begins to be 
submerged and is trapped into the material that is pil- 
ing up onto the neutron star surface. The matter outside 
the accretion shock is still falling with constant accretion 
rate, but it is the fiuid inside the envelope, with its much 
lower downward velocity, which is nevertheless responsi- 
ble for the magnetic field submergence. After 100 ms, 
the submergence is completed, with a maximum mag- 
netic field strength :^ 5 x 10^^ G. At t = the magnetic 
loop was immersed in a medium of density 10^ g cm~^, 
the infalling matter just behind the reverse shock had a 
density ^ 10^ g cm~^ but after 100 ms the magnetic field 



has been submerged into matter at density 



5 X 10^ g 



For comparison, in Figure [4] we show the results of a 
similar simulation but within the QHEE scenario. At 
t = the magnetic field loop is immersed within the 
QHEE. The initial, strongly convective phase seen in the 
free-fall scenario is absent. Some weak convection is nev- 
ertheless present in the upper part of the envelope, whose 
height initially oscillates {t = 1 panel in Figure [S. Due 
to the smoother evolution, compared to the free-^11 sce- 
nario, even at this early time the magnetic loop has al- 
ready been compressed by the infalling envelope but is 
not strongly disrupted as in the free-fall scenario. After 
100 ms, the QHEE has recovered its initial shape and 
height, with an equilibrium shock radius rgh — 4 x 10^ 
cm, and the magnetic field is seen to be submerged. Al- 
though the accretion velocity within the envelope is much 
lower than above the shock, it is enough to completely 
submerge the magnetic loop. 

In the free-fall scenario the initial smashing of the mag- 
netic loop by the reverse shock is immediately counter- 
acted by the bouncing of the matter off the neutron star 
surface which drags the field with it, upward, within the 
forming envelope. The field submergence in the QHEE 
scenario at early times shows that it is the slow settling of 
the envelope material, under continuous accretion, which 
is actually responsible for the field submergence. In the 
free-fall scenario the field submergence follows the same 
route as in the QHEE scenario once the envelope stabi- 
lizes and its matter slowly settles onto the neutron star 
surface. The final geometry of the magnetic field is seen 
to be very similar under both free-fall and QHEE initial 
conditions. 

For the lower accretion rates, 10 and Irho, the evolu- 
tion is very similar but the accretion shock reaches much 
larger heights. We choose to illustrate only the 100 mo, 
in Figures [3] and [4] since the lower rate cases involve much 
higher columns. We show in Figure |5] the radial profiles 
of the magnetic field for our three highly hypercritical 
accretion rates, in the free-fall scenario, corresponding 
to the times at which the quasi-hydrostatic equilibrium 
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Fig. 3. — Density maps (color) with magnetic field contours superimposed, at three instants, t = 0, 1, and 100 ms, for a 2D MHD 
simulation in the free-fall scenario, with rh = 100 rho- The imposed level of refinement was 4, with 3 blocks in x-direction and 39 in 
the y-direction, which results in 192 x 2496 effective zones in the computational domain. Note how the shock rises progressively and the 
magnetic field is submerged close to the neutron star surface. The anchor points of the magnetic loop are still visible in the final snapshot. 
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Fig. 4. — The same as in Figure [3] at the same times t = 0, 1, and 100 ms, but for the initial atmosphere in quasi-hydrostatic equilibrium 
(QHEE). After the initial transient, the fiow settles also to a quasi-static atmosphere, with the magnetic field being strongly confined to 
the vicinity of the neutron star surface. The anchor points of the magnetic loop are also visible in the final snapshot. 
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Fig. 5. — Radial profiles of the magnetic field strength, once 
the stationary state has been reached, in the free-fall scenario, for 
777, = 777,Q at time t = 600 ms, m = 10 mo at time t = 300 ms, and 
777, = 100 mo at time t = 100 ms. The horizontal averages of the 
field strength as a function of the height above the neutron star 
surface are plotted for each case (labeled). 



has been reached. In the three cases the submergence of 
the magnetic field is clearly seen. 

It is thus clear from the above results that for strongly 
hypercritical mass accretion rates, both approaches, the 
free fall and QHEE initial conditions, lead to the same 
conclusion, namely, that the field is promptly submerged 
by the sheer force of the settling fiuid in the envelope. 
Thus we believe this to be a robust result refiecting the 
actual fate likely to befall an initial magnetic field an- 
chored on the neutron star once a quasi-stationary set- 
tling envelope has been established. 

3.1.1. 3D simulation 

In all the high accretion rates described above we ob- 
serve the submergence of the magnetic field. These were, 
however, 2D geometries and, since the magnetic field be- 
havior is an intrinsically 3D phenomenon, one may won- 
der how much the addition of the third dimension degree 
of freedom would alter the results. Moreover, in a 3D 
environment the accreting matter has the possibility to 
slip around the magnetic field loop, hence potentially re- 
ducing the pressure it exerts on it. To tackle this issue, 
we performed a 3D MHD simulation for the free-fall sce- 
nario. From the available computing power we were lim- 
ited to simulate only a short column, hence we consider 
only our highest accretion rate, 100 mo- 

In Figure |6] we show four snapshots of the time evo- 
lution of the magnetic field iso-surfaces until the quasi- 
statonary state is reached. Notice that the fiuid reaches 
the stationary state in 60 ms, i.e. in a shorter time than 
in the corresponding 2D simulation. Overall, the evolu- 
tion is very similar to the corresponding 2D simulations. 
After 1 ms the magnetic field loop has been smashed 
by the accretion shock and the field is beginning to be 
spread by the violent fiuid motion. The accretion shock 
rapidly progresses upward and reaches its final height, 
rsh — 4 X 10^ cm as in the 2D case, in about 20 ms. No- 



tice, at that time, that the magnetic field is also being 
dragged upward by the fiuid motion but is not able to 
reach the location of the shock. In no moment will the 
field be able to rise above the forming envelope. However, 
when the QHEE has built up and the initial strongly con- 
vective regime has died away the magnetic field has been 
forced downward by the slow settling of matter in the 
envelope and is submerged inside the denser layers onto 
the neutron star surface. Notice that the magnetic field 
has been spread all over the whole bottom area of the 
simulation box. 

In Figure [7| we show the radial profiles of the density, 
pressure, velocity and magnetic field in the last stage. 
The density and pressure profiles are essentially identi- 
cal to the 2D ones shown in Figure |2] with the clear piling 
up of matter onto the neutron star surface. The velocity 
profile, however, is much less noisy in 3D than in 2D. 
These three, p, P, and profiles are also in close agree- 
ment with the analytical solution. The magnetic field 
strength profile clearly exhibits the submergence within 
the high density region. We thus conclude that at this 
level of mass accretion, no significant differences in mor- 
phology are apparent by considering the 3D nature of 
the problem, and our previous statements concerning the 
submergence of the field from the 2D calculations remain 
valid. 

As a further energy check, the neutrino luminosity in- 
tegrated over the whole domain, at t = 60 ms when 
the quasi- st at ionnary state has been reached, is Lj, :^ 
2.3 X 10^^ erg s~^ but the emission is strongly concen- 
trated within the first two kms above the stellar surface 
due to the strong temperature dependence of the emissiv- 
ity. The gravitational energy liberated by the accretion 
is ~ 0.3 X mc^ 2.1 X 10^^ erg s~^ and is hence emitted 
essentially in its entirety in neutrinos. 

Once the quasi-stationary state has been achieved, at 
t = 60 ms, the adiabatic indices connected with the 
sound velocity and with the system energy are: 7c = 
1.35 and 7e = 1.34. Then the adiabatic and radia- 
tive gradients are: Vad = 1 — I/Tc — 0.26, Vrad = 
{d\nT/d\nP) :^ 0.25. The system is, hence, convec- 
tively stable, but only marginally. 

3.2. Weakly hypercritical accretion rates 

We can now investigate what happens with the 
magnetic loop anchored onto the neutron star sur- 
face for lower accretion rates. To answer this, 
we carry out 2D simulations for the cases m = 
(10~^, 10~^, 10~^, 10~^) • mo- For these lower accre- 
tion rates it is computationally impossible to follow the 
expansion of the accretion shock because it reaches to ex- 
tremely large heights, leaving the computational domain, 
and, moreover, the relaxation time of the envelope be- 
comes exceedingly long. Attempting to follow the evolu- 
tion of a free-fall scenario by allowing the accretion shock 
to leave the integration domain would be physically, and 
numerically, inconsistent with our setting: it not possible 
to impose an upper boundary condition which incorpo- 
rates in any reasonable manner the strongly turbulent 
motion of the initial phase of envelope formation. Due 
to this, we are forced to restrict ourselves to the QHEE 
scenarios. The results of o ur s tudy of the highly hyper- 
critical accretion rates in ^13. 11 showed that the field sub- 



mergence is occurring in a very similar manner in both 
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Fig. 6. — Iso-surfaces of magnetic field at times t=0, 1, 20, and 60 ms, for the 3D simulation in the free-fall scenario, and an accretion 
rate m = 100 mo- We show color maps of density and pressure in the xy and yz faces, respectively. In the two upper panels only the lO-*^-*^ G 
iso-surfaces are visible, the magnetic field being, however, stronger within these surfaces. In the two lower panels the 10-*^ -"^ G iso-surface is 
rendered with high transparency allowing to clearly see much stronger magnetic field regions. The end result of magnetic field submergence 
due to the strong accrretion is the same as in the 2D simulations. 
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Fig. 7. — Radial profiles of averaged density, pressure, velocity 
and magnetic field for the 3D model depicted in Figure [6] once the 
stationary state has been reached, at t = 60 ms. Comparison with 
Figure [2] shows that the solution in terms of position of the shock 
front and values is the same as in 2D. 

types of scenarios and that, in the long term, it is manly 
due to the settling of matter onto the neutron star sur- 
face. This gives us confidence that limiting ourselves to 
the QHEE scenario for the weakly hypercritical rates will 
give us representative results. 

We simulated an accretion column with height /S.y = 
4x 10^ cm, using the upper boundary condition with mat- 
ter injection at the same m and the corresponding p and 
V as given by the analytical solution for a QHEE. The ini- 
tial magnetic field configuration was that of a magnetic 
loop, as described in §|2j immersed within a QHEE enve- 
lope Each simulation exhibited a short initial relaxation 
phase but the system rapidly reached stationary density 
and pressure profiles. In Figure [8] we show color maps 
of density, with magnetic field contours superimposed, 
for these four accretion rates at times where a stationary 
state was well established in the density profile. At the 
highest accretion rate, rh = 0.1 mo, in 100 ms the mag- 
netic loops is clearly submerged within the high density 
bottom of the envelope. This result is similar t o wh at we 
obtained at the highly hypercritical rates in § |3.1[ 

At lower rates, m <C 0.1 rho, however, the situation 
becomes more complex. As seen in Figure [Sj the mag- 
netic field is mostly confined in the high density region 
but is able to extend to slightly larger heights than be- 
fore. The evolution of our lowest accretion rate scenario, 
m = 10~^ mo, is depicted in greater detail in Figure [oj a 
simulation we could pursue for 400 ms (after that time 
the turbulence reached the upper boundary of our in- 
tegration domain and confiicted with the infiow condi- 
tion, so the simulation was stopped). The left panels 
show that, after an initial strong compression, the field 
is slowly rising: at 100 ms it is essentially confined within 
the high density region but, as matter piles up onto the 
neutron star surface and this high density region expands 
upward, the magnetic field follows it. The right column 



in Figure [9] illustrates the pattern of convective cells with 
very high velocities in the upper low density regions that 
contrast with the much slower motion within the much 
denser matter accumulating onto the neutron star sur- 
face. The size of the convective cells is expected to be of 
the order of the pressure scale- height. Hp = —dr/d\ogP^ 
which is ~ 10 km in the envelope but only ~ 1 km in 
the forming crust. These two distinct patterns, and their 
length-scales, are clearly seen in Figure [Toj where the ab- 
solute value of the vertical velocity is plotted as a func- 
tion of radius, and are conserved during the evolution in 
spite of a more than one order of magnitude decrease in 
velocity. It is this clear difference in the two convection 
patterns, large, low density cells with high velocites on 
top of small, high density cells with low velocities that 
keeps the crushed magnetic field confined into the high 
density region of the forming new crust. 

4. DISCUSSION AND INTERPRETATION 

We have continued our study of hypercritical accretion 
onto newborn neutron stars, in particular with the aim 
of studying the possible submergence of the magnetic 
field for various accretion rates. This is particularly rele- 
vant in the context of making the neutron star eventually 
invisible as a pulsar following the supernova explosion. 
Extending our exploration of parameter space to wide 
columns spanning a significant fraction of the stellar sur- 
face, we have paid special attention to the formation of 
the quasi-steady atmosphere that forms after the reverse 
shock has reached the hard surface of the star, and how it 
affects an initial loop of field anchored in it. We find that 
whether an initial condition in free fall, or an analytical 
atmosphere in quasi-hydrostatic equilibrium, is initially 
used, the field is rapidly submerged into the forming 
high-density new crust for strongly hypercritical rates, 
above our fiducial value of M = Mq = 350 Mq yr~^. 
The comparison between results for different initial con- 
ditions is important because it is extremely challenging 
to follow the evolution of the system, in particular the 
shock front, for lower accretion rates. This is due to the 
scaling of the equilibrium position of the shock, given in 
Equation |8) that shows the shock expands to very large 
heights at small accretion rates. We have further car- 
ried out a calculation in three dimensions, to investigate 
whether the fiow exhibits qualitative changes, in partic- 
ular regarding the submergence of the field, as compared 
to the 2D simulations. Finding none, we conclude that 
2D studies are able to correctly capture the main features 
of the field's behavior in this respect. 

As the accretion rate is lowered, reaching down to 
10~^ Mo ~ O.O35M0 yr~^, the field is progressively less 
affected by the infall, and the initial loop reaches upwards 
of the neutron star surface after the initial transient, to a 
height of a few kilometers. Intringuingly, the pattern of 
convective cells helps explain why the field is not dragged 
to greater heights. Close to the surface of the star, a high- 
density region where matter accumulates is present and 
this is also where most of th e magnetic field remains con- 
fined. As described in §3.2[ there is a sharp difference in 
the convection patterns in the high density region, with 
cell sizes ~ 1 km, compared to the overlaying envelope, 
with cell sizes ~ 10 km. Thus, the envelope plasma cir- 
culations are unable to reach into the region of high field 
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Fig. 8. — Color maps of the density with magnetic field intensity contours superimposed, for the four weakly hypercritical accretion rates 
in the initially QHEE scenario: 10~ and 10~ mo at time t = 100 ms, 10~^ mo at time t = 200 ms, and 10~^ mo at time t = 400 ms. The 
submergence of the field is now clearly affected by the assumed accretion rate. 
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Fig. 9. — Evolution of the accretion fiow in the QHEE initial condition scenario with m = 10~^ mo at 6 different times, as labeled. Left 
panels: color maps of the density with magnetic field intensity contours superimposed. Right panels: color maps of the velocity with arrows 
showing its direction. After an initialy transient in which the field is rapidly compressed, it rises to a height of 2^ 4 km (radius ~ 14 km). 
The change in the pattern of the fiow across the boundary from low (at top) to high (bottom) density is apparent. 
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Fig. 10. — Profiles of the absolute value of the vertical velocity, 
\vy\^ in the rn = lO^'^mo scenario at t = 10 ms (dotted curve) 
and t = 400 ms (solid curve). The position of the shock from the 
infalling matter at 10 ms is labeled. The regular fiuctuations in 
velocity are due to the convective motions in the fiow. The size of 
the cells can be readily measured and the widely different lengh- 
scales between the high density forming crust and the envelope are 
clearly seen. 



and drag it to larger altitudes. Notice that even at the 
lowest accretion rate we considered, the matter pressure 
is much larger than the magnetic pressure and the field 
is continuously dragged by the motion of the plasma. 

Our distinction between highly and weakly hypercrit- 
ical accretion rates is purely due to our computational 
limitations. However, we do find a physical distinction 
in the field submergence that is simply due to the evo- 
lution of the convection in the envelope and the forming 
crust. In the highly hypercritical cases of §3.1| and the 
highest rate, O.lrho, of ^3.2 our simulation lasted long 
enough for convection in the envelope to essentially disa- 
pear, while at the lower rates significant convection was 
always present. In the former cases the only remaining 
fiuid motion is the slow settling of the accreting mat- 
ter onto the neutron star and, since the magnetic field is 
frozen into the plasma, the field is completely submerged 
into the high density region of the forming crust. This 
result is still clearly seen in the O.lmo case of Figure [8) 
In contradistinction, in the latter cases continued con- 
vection maintains mixing of the magnetic field, but the 
different patterns of convection in the forming crust com- 
pared to the envelope traps most of the magnetic field in 
the high density regions and the field present in the en- 
velope is strongly reduced compared to its initial value. 
The evolution displayed in the left panels of Figure [9] 
clearly shows the slow submergence of the field in the 
high density region as matter piles up onto the neutron 
star and the continuous presence of a residual field in the 
bottom of the envelope in the interface region between 
the two regions of different convection patterns. Notice 
that the apparent rise of the magnetic field in this Fig- 
ure [9] is actually just the rise of the height of the building 
new crust while the residual field never extends into the 
envelope much beyond the interface region. If this accre- 
tion phase last for several hours, most of the magnetic 
field should eventually be pushed deep into the crust. 

The critical factor separating a total versus a partial 



submergence of the initial magnetic field hence seems to 
be the time scale for convection in the accretion envelope 
to disappear. Once convection stops the magnetic field is 
unavoidably submerged into the neutron star curst. As 
described by Fryer et al.j (^1996), at low accretion rate 
the development of a stationary non-convective envelope 
may take more time than the duration of the hyper- 
critical accretion phase. Supernovae in which a strong 
reverse shock was produced and in which the accreting 
envelope reached a convectively stable profile are likely 
to produce a neutron star with a vanishingly small sur- 
face magnetic field. In contradistinction, supernovae that 
produced a weak reverse shock whose accreting envelope 
remained convective can be expected to produce weakly 
magnetized neutron stars. In the case where no reverse 
shock developed, the neutron star will simply maintain 
the magnetic field it had at its birth. 

5. CONCLUSIONS 

We have extended our previous study of Paper I of 
the effect of post-supernova hypercritical accretion on 
the magnetic field of the new-born neutron star. Our 
study is still restricted to cases where the accreting mat- 
ter pressure dominates over the magnetic pressure, so 
that the magnetic field is only a "victim" of the ac- 
cretion. We found that the critical factor determining 
the fate of the magnetic field is whether the develop- 
ing quasi-hydrostatic equilibrium accreting envelope will 
reach of state of convective stability or not. Once convec- 
tion stops the magnetic field is unavoidably pushed into 
the forming new crust and a non magnetized neutron 
star will result. If the acreting envelope remains con- 
vective, most of the initial magnetic field is still pushed 
into the forming crust but the fiuid motions, due to the 
diferent patterns of convection in the crust and the sur- 
rounding envelope, maintain a residual field protruding 
out of the star. In this latter case, the resulting exter- 
nal field certainly has a very complicated structure and 
anisotropy in the accretion, as, e.g., from rotation, will 
also increase the compexity of the final field. The dipo- 
lar component of the final field may be several orders 
of magnitude smaller than the initial one but localized 
spots with a much stronger field seem unavoidable. 

Non-convective envelopes are rapidly obtained only at 
extremely high accretion rates, larger that ~ IOMq yr~^ 
in our study. It is not clear how many new-born neu- 
tron stars undergo such events, and how many of them 
would survi ve without collapsing into a black- hole (Fryer 
et al.| 1996). Our conclusion is that there may be only 
very few young neutron stars that have been totally de- 
magnetized by hypercritical accretion and they are likely 
very massive ones. 

At lower accretion rates, field submergence is only par- 
tial and the fraction of new-born neutron stars having 



passed throu gh such events is likely significant, p^ryer 
et al. (1996) were very careful in the treatment of neu- 
trino heating and found that when M > 0.1 Mq yr~^ a 
neutrino driven explosion may result, expelling the en- 
velope while at lower rates steady accretion is achieved 
with, eventually, a non-convective envelope. As the ac- 
cretion rate decreases one may reach a point where the 
residual magnetic field may begin to play a dynamical 
role. The final outcome of such evolution is impossible 
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to treat numerically and can only be guessed at now, but 
a weakly magnetized neutron star, with a complicated 
field geometry seems very likely. 

The CCOs are the primary candidates for neutron stars 
having undergone such a field submergence and they con- 
stitute lialf_of_th£^n£w^ of young neutron 



stars: 



Popov & Turolla (2012) list a dozen of pulsars 



with age < lU kyrs, three ot them b eing CCOs, i.e. a to 



tal of nine non-CCOs pulsars, while Gotthelf & Halpern 
(|2008) present eight CCOs, all having ages below 10 
kyrs to which one may add the yet undiscovered neu- 
tron star in the remnant of SN 1987A. The three CCOs 
with measurements of spin periods P, and also deriva- 
tives P from which surface dipolar field strengths Bs can 
be inferred, are PSR J1852+0040 in Kesteven 79, with 
a period P = 105 ms and = 3.1 x 10^° G (lHalpern| 
etar][2QQ7| , and PSR J1210-5226 in PKS 1209-51/52, 
with P = 424 ms and Bs = 9.8 x 10^° G and PSR J0821- 
4300 in Puppis A with P = 112 ms and Bs = 2.9 x 10^° G 
(Gotthelf & Halpern to be published in ApJ). 

The unavoidable prediction of a submerged magnetic 
field is the subsequent growth of th e surface field by back - 
diffusion of the internal field (Muslimov fc Page 1995). 
Hence, burial of the field in a newborn neutron star may 
manifest itself as a delay in the star turning on as a classi- 
cal radio pulsar. Possibl e evolution scenarios were stud- 
ied in ID models in Geppert et al. (1999) and |Ho| (12011 ) 



and recently by Vigano Pons (2012) in a 2D scenario. 
A growing magnetic held also natura lly manifests itself 



through a b raking index inferior to 3 (Muslimov & Page 
1996 1999) and may lead to peculiar evolutionary prop- 



erties, such as those a lready observed in PSR J 1734-3333 
dEspinoza eFar||2011 ). 

in the present work we were still dealing with mag- 
netic field strengths, and accretion rates, in which the 
magnetic forces are negligible, even though the FLASH 
code does take them into account. More interesting ef- 
fects can be expected in the case of a magnet ar size field 
( Thompson fc Murray||2001 ) and will be considered in a 
future paper. 

Finally, among the many simplifications of our model, 
the assumption of non magnetized accreted matter is 
clearly disputable. In case the neutron star magnetic 
field is of fossile origin, the accreting flow is likely bring- 
ing in some significant magnetic field that may just re- 
place the submerged one. One may conceive that accre- 
tion will hence complicate the geometry of the surface 
field without necessarily reducing its strength, but pos- 
sibly reducing the strength of the dipolar component. 
However, turbulence in t he accretion envelope can act as 



a diamagnetic medium ( Vamshtem fc Zerdovich] 1972) 
which may also result in a strong suppression of the sur- 
face magnetic field. Such an evolution certainly deserves 
further study. 
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